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Abstract 



We report results from shell model studies of muon capture on '^Ca to 
low-lying levels of ^'^K. We discuss the comparison between calculated capture 
rates, measured capture rates and analogous transitions in (e,e') scattering 
in terms of the particle-hole structure of the ^'^Ca-'^'^K nuclei. We highlight 
the ^''Ca(0"'", 0)— >^'^K(0~, 2626) axial charge transition and its sensitivity to 
the induced pseudoscalar coupling gp of the proton's weak interaction. In ad- 
dition, we address the hindrance of unique first-forbidden transitions due to 
particle-hole interactions and the emergence of allowed Gamow- Teller transi- 
tions due to ground state correlations. Lastly, we examine the longitudinal 
alignment of '^''K recoils following muon capture, and discuss this possibility 

for independently determining the induced coupling gp. 
23.40.-S, 23.40.Hc, 27.40.+Z 
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I. INTRODUCTION 



Studies of muon capture on complex nuclei have spanned nearly fifty years. The early 
work [1-3] was devoted to establishing the universal V-A character of nuclear muon capture, 
while later work has focused on induced currents [4], second-class currents [5,6] and non 
V-A interactions [7]. Muon capture is also a valuable window on particle- hole excitations 
and spin-isospin modes in nuclei [8] , it complementing the information that is obtained from 
bet a- decay, electron scattering and charge exchange reactions. 

Of particular interest in nuclear muon capture is the induced pseudoscalar coupling Qp 
of the proton's weak interaction. For the free proton, the couphng gp is predicted to few- 
percent accuracy by symmetry arguments and thereby represents an important test of low 
energy QCD [9-11]. For the bound proton, its medium modification is sensitive to effects 
that range from pion exchange currents and A-hole excitations to partial restoration of chiral 
symmetry [12-14]. Indeed, some suggestions of a large A-dependent renormalization of the 
induced coupling gp have been published in the literature [12,15]. 

Measday and Stocki [16] have recently published new experimental results on 7-ray spec- 
tra from muon capture on ^°Ca that include determinations of partial rates and rate limits 
for numerous ^'^Ca(/x,i/)^°K transitions. Their data are noteworthy as the heaviest nucleus on 
which (/i,z^) reactions have been clearly identified. Their data show evidence of excitations 
that range from allowed Gamow- Teller transitions to highly ^-forbidden transitions. 

Most important is their observation of the ^°Ca(0+,0) — > ^°K(0~,2626) axial charge 
transition. In nuclear beta decay such 0"'"<<-^0~ transitions have been extensively studied - 
both experimentally and theoretically - and nowadays represent our clearest evidence for 
exchange currents in complex nuclei. In muon capture the only known example of a 0"'"<->0~ 
first forbidden transitions was the celebrated ^^0(0"*", 0)— i^^^N(0~, 120) transition, it having 
been studied by many authors in the context of the coupling gp. In light of suggestions of an 
A-dependent renormahzation of the induced coupling gp, the ^'^Ca(0+,0) — > ^°K(0~,2626) 
transition represents a valuable data point for experimentally accessing gp in medium-weight 
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nuclei. 

Also interesting are two unique first forbidden transitions and one allowed Gamow Teller 
transition: 40Ca(0+,0) ^ ^°K(2-,800), ^°Ca(0+,0) ^ ^0^(2", 2047) and ^0Ca(0+, 0) ^ 
^°K(1+, 2290). One motivation for investigating unique first-forbidden transitions stems 
from the large quenching of such transitions in nuclear /3-decay, such quenching arising from 
the coherent effects of the repulsive T=l particle- hole interaction. One motivation for inves- 
tigating allowed Gamow- Teller transitions stems from their absence in a simple closed-shell 
description of doubly-magic ^'^Ca, their emergence arising from the multi-particle, multi- 
hole admixtures in the ^°Ca ground state. In addition, the '^^Ca{0^ ,0) — > **°K(2~,800) and 
^°Ca(0+, 0) — > ^°K(2~, 2290) transitions have well known analogs in (e,e') inelastic scattering 
[17-19] and charge exchange reactions [20,21]. The comparison between u) capture rates 
and analogous (e,e') transitions thus provides a valuable cross-check for the ^°Ca(/i,i/)'^°K 
data and a valuable comparison for the ^°Ca(//,i/)^'^K calculations. 

The ^°Ca(//,i/)^°K data are also interesting in the broader context of the shell structure 
of the A 40 region. Warburton and co-workers have published a series of articles [22-26] 
on mass 34 < A < 50 nuclei that include comprehensive studies of AJ^ — 0~ P decays and 
their enhancement due to exchange currents effects, A J'^ — 2~ P decays and their hindrance 
due to core polarization effects, and the evolving nuclear structure of the A ~ 40 mass 
region. Weak transition rates on doubly magic ^°Ca are the archetypal benchmark for such 
model calculations. 

Herein we report results from shell model studies of muon capture on "^°Ca. In Sees. II 
and III respectively we describe the details of the shell model calculations and the capture 
rate calculations. In Sec. IV we compare the calculated (//,i^) transition rates with the recent 
data of Measday and Stocki [16] and the earlier data of Igo-Kemenes et al. [27] in the context 
of the particle-hole structure of the ^°Ca-^°K nuclear wavef unctions. In particular, we discuss 
the emergence of the (0+, 0)^(1+, 2290) allowed Gamow-Teller transition, the hindrance of 
the O"*"— >-2~ unique first-forbidden transitions, and the sensitivity to the induced pseudoscalar 
couphng Qp of the ^°Ca(0"'", 0)^''°K(0~, 2626) axial charge transition. Finally, in Sec. V we 
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briefly consider the longitudinal alignment of °K recoils following ^°Ca capture, and the 
dependence of the alignment on the coupling gp. 

II. SHELL MODEL CALCULATIONS 

The simplest desription of muon capture from the ^'^Ca ground state to the low-lying 
states involves a fully filled sd-shell initial state and various (sd)"^(pf)^ particle- hole 
final states. This model infers a significant number of forbidden transitions but no allowed 
transitions as all spin-orbit partners in the ^'^Ca ground state are fully filled. It implies 
that data for allowed transitions will have interesting sensitivity to the multi-particle, multi- 
hole admixtures in the "^"Ca g.s. wavefunction and data for forbidden transitions will have 
interesting sensitivity to the particle-hole interaction in the A ~ 40 mass region. 

The ^'^K level diagram of Fig. 1 reveals one negative parity quadruplet with spin-parities 
J'^ = 4~,3~,2~ and 5~ and excitation energies 0,30,800 and 891 keV and another neg- 
ative parity quadruplet with spin-parities J'^ — 2~,3~,1~ and 0~ and excitation energies 
2047, 2070, 2104 and 2626 keV. The quadruplets are consistent with low-lying negative parity 
states involving ((i3/2)~^(/7/2)^ particle-hole configurations, which yield the J'" = 2~- 
5~ quadruplet, and (^3/2)^^(^3/2)^ particle-hole configurations, which yield the J'^ — 0~- 
3~ quadruplet. In addition, the four lowest-lying positive parity states, with spin-parities 
J'^ = 0+, 2+, 3+ and 1+ and excitation energies 1632, 1959, 2260 and 2290 kcV, were identified 
by Davis et al. [28] as consistent with expectations of low- lying {ds/2)~'^{f7/2)^ excitations 
with the two /7/2-particles coupled to J = and the two c?3/2-holes coupled to J = 0- 
3. At higher excitation energies the increasing level density makes meaningful assignments 
between model configurations and experimental levels much more difficult. 

Below we describe the model spaces and residual interactions we have employed in our 
calculations of the muon capture rates and the ''^K recoil alignments. Our basic philosophy 
was to utihze well-known, well-tested models "off-the-shelf", i.e. not fine tuning any model 
parameters. Our model calculations were performed utilizing the OXBASH code [29] . 
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A. ^ = 40 model spaces and residual interactions 

Our primary calculations of ^°Ca(//, i/)^'^K capture were conducted using the sd-pf model 
space and the WBMB interaction [30]. The WBMB interaction comprises the sd interac- 
tion of Wildenthal [31], the pf interaction of McGory [32], and a modified Millener-Kurath 
interaction for the cross-shell matrix elements [33]. As discussed in Ref. [22], the sd-pf shell 
energy gap was empirically adjusted to reproduce the A = 35-41 binding energies and the 
(^3/2 /7/2'*-^ 0^3/2/7/2, ds/2P3/2'^ds/2P3/2 two-body matrix elements were empirically adjusted to 
reproduce the low-lying A — AO particle-hole states. 

Two versions of WBMB calculations were carried out. In both versions of the WBMB 
model the negative-parity states comprised one particle in the fp shell and one hole in 
the sd shell, i.e. {sdy^{pfy configurations, and the positive-parity states comprised two 
particles in the fp shell and two holes in the sd shell, i.e. {sdY'^{pfY configurations. However, 
to enable the investigation of ground state correlations we employed ^°Ca wavefunctions 
with both a pure {sd)"^ configuration, our so-called hu WBMB calculation, and a mixed 
{sdY'^+{sd)'^'^{pfY configuration, our so-called (0+2) hto WBMB calculation. We caution 
the reader the WBMB interaction was originally developed and generally applied in pure 
N hw model spaces, although the effects of (0-1-2) hw mixing on hypothetical ^°K— >^°Ca 
beta-decays was discussed in Ref. [23].^ 

We also performed calculations that employed a smaller 0^3/2-/7/2 model space but allowed 
a richer 0-8 hw particle-hole spectrum. These calculations used the single particle energies 
and two-body matrix elements of Sakakura, Arima and Sebe [34] which is hereafter denoted 
as the SAS interaction. This interaction was used in Ref. [34] to study the effects of multi- 
particle, multi-hole admixtures in 0"^ levels of ^'^Ca. While the smaller (^3/2-/7/2 space is 
obviously unable to fully describe all low-lying "^^K particle-hole levels - it missing the fzh- 



^The correct treatment of (0 -|- 2)%u mixing is a well-known difficulty in shell model calculations 
See Ref. [22] for a discussion of (0 + 2)hu} mixing in the A ~ 40 mass region. 
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P3/2 states - it renders some insight into the multi-particle, multi-hole components in the 
^°Ca ground state wavefunction and their influence on the (/i, u) transitions. 

B. excitation energies and wavefunction configurations 

The WBMB calculations exactly reproduces the excitation energies of the lower-lying 

= 4~, 3^, 2^, 5^ quadruplet and the ^'^K higher-lying J"^ = 2~, 3^, 1^, 0^ quadruplet due 
to the aforementioned adjustments of the (^3/2/7/2^ 0^3/2/7/2, d^/2P^/2'^d^/2Pz/2 two-body 
matrix elements. Within the WBMB model the J'^ = 4~, 3~, 2~, 5~ quadruplet has d~J^-fy^ 
occupancies that range from 92 to 98% and the J'^ = 2~, 3~, 1~, 0~ quadruplet has d'^j2~P\/2 
occupancies that range from 75 to 99% (the occupancies are consistent with the compilation 
of spectroscopic factors in Endt et al. [35]). By comparison the WBMB model yields low- 
lying positive parity states that are typically 1 MeV under-bound and have somewhat 
exaggerated energy spacings ~ although the model does reproduce the J"" = 0"*", 2+, 3"*", 1"*" 
level ordering. While the J'^ = 0+, 2+, 3+ states have dominant d'^J^f 7/2 character the J'^ — 
1"*" state is a more fragmented mixture of the permissible 2p-2h configurations. 

The SAS model reproduces the excitation energies of the lower-lying negative-parity "^^K 
quadruplet with spin-parities J'^ — 4~, 3~, 2~, 5~ to accuracies of roughly 200 keV. The cor- 
responding SAS wavefunctions have their largest contributions from d^J^^f 7/2 configurations 
(43-55%) but with substantial contributions from d'^j^f 7/2 configurations. (31-40%). The 

— 0+, 2+, 3"^, 1+ level ordering is reproduced by the SAS model although the J'^ — 0+, 2+ 
states are over-bound by roughly 1 MeV. All four positive parity states have their largest 
contributions from d'^J^fy^ configurations (42-61%) but with substantial contributions from 
^3/2/7/2 configurations (34-40%). Of course the SAS model and (^3/2-/7/2 space is unable 
to account for the higher-lying negative-parity quadruplet with its dominant (^3/2-^3/2 
character. 
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C. Ca ground state wavefunction and np-nh admixtures 

The various calculations are quite different in their descriptions of the ^°Ca g.s. While 
the hcu WBMB calculation assumes a simple closed-shell ^°Ca g.s. the closed-shell occu- 
pancies are 72% in the (0-1-2) hu WBMB calculation and 46% in the SAS calculation. In 
addition, the (0-1-2) ?iw WBMB calculation and SAS calculation have qualitatively different 
np-nh wavefunction admixtures. For example, in the SAS calculation the leading np-nh 
component is the (^3/2-/7/2 configuration with a 35% occupancy whereas in the (0-1-2) Jlu 
WBMB calculation the large np-nh fragmentation yields a d^J^-fy^ configuration with only 
2% occupancy. 

The different descriptions of the ^°Ca g.s. - in particular the differences in closed-shell 
occupancies and multi-particle, multi-hole admixtures - is the leading cause of the model- 
to-model variation in the calculated capture rates. It is worthwhile noting the negative- 
parity transitions and positive-parity transitions are endowed with qualitatively different 
sensitivities to the np-nh admixtures in the ^°Ca ground state. This difference originates as 
negative-parity transitions are dominantly one-body transitions connecting the ^°Ca closed- 
shell wavefunction component and the Ip-lh wavefunctions while the positive-parity 
transitions are dominantly one-body transitions connecting the ^^Ca 2p-2h wavefunction 
components and the ^'^K 2p-2h wavefunction. Consequently, the negative-parity transitions 
are mostly sensitive to the magnitude of the ^°Ca 2p-2h admixtures while the positive-parity 
transitions are additionally sensitive to the configurations in the ^'^Ca 2p-2h admixtures. 

III. MUON CAPTURE CALCULATIONS 

In calculating the muon capture observables and nuclear matrix elements we used the 
formalism and notation of Walecka [37] and Donnelly and Haxton [38] (for details see Ref. 
[39]). Note that in ^°Ca capture each 0"^^JJ transition involves a single spin-parity mul- 
tipole and all observables are governed by the nuclear matrix elements of two electroweak 
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operators, i.e. Cj — Aij and 7^'_7j""f principle these nuclear matrix elements can 
involve the contributions of one-body, two-body and many-body weak nucleonic currents. 
Herein we assume these operators can be represented by the sum of the A one-body weak 
nucleonic currents, i. e. the impulse approximation. In such circumstances the nuclear matrix 
elements oltj — Mj and 'fei_q-p<^9 expressible as sums of products of single-particle 
matrix elements (SPMEs) and one-body transition densities (OBTDs). 

In calculating the nuclear matrix elements we assumed harmonic oscillator wavefunc- 
tions with an oscillator parameter h = 1.94 fm. In addition, we assumed a constant muon 
wavefunction over the nuclear volume with |0is|o„ = R|0is(O)p where 0is(O) is the muon 
wavefunction for a point-hke nucleus and R=0.44 is the reduction factor for the ^'^Ca nucleus 
[37]. In calculating the muon capture rates we generally fixed the numerically values of the 
weak vector, weak magnetic, weak axial and induced pseudoscalar couplings constants at 

— 1.0, Qm — 3.706, Qa — 1.27 and Qp — 8.23. While the coupling Qp is defined at = -|-0.9m^, 
the other couplings are defined at = and therefore were scaled assuming a dipole form 
factor with = 0.73 GeV^. The momentum transfer q was computed via q + q^/2Mt = rUf^ 

— AE — eb where is the target mass, m^ is the muon mass, is the /i" binding energy, 
and AE is the energy difference of the ^'^Ca-^^K nuclear states. 

As mentioned above we employed harmonic oscillator wavefunctions in calculating nu- 
clear matrix elements - such wavefunctions poorly reproducing the true asymptotic behavior 
of nuclear radial wavefunctions. We note the effect of different choices of radial wavefunc- 
tions in first forbidden beta-decay in A ^ 40 nuclei was studied by Warburton et al. [23]. 
These authors reported only small differences between the relevant nuclear matrix elements 
obtained with harmonic oscillator wavefunctions and Woods-Saxon wavefunctions. In par- 
ticular, the d3/2<->^P3/2 matrix elements relevant to (0"*", 0)— >(0~, 2626) capture were found 
to differ by roughly 5-10% (of course it is possible that their results are not apphcable at 
the higher momentum transfer of the muon capture process). In addition, the sensitivity of 
weak matrix elements and muon capture rates to different radial wavefunctions was studied 
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in A — 12-32 nuclei by Kortelainen et al [36]. Again, in most cases these authors observed 
only minor differences between results derived with harmonic oscillator wavefunctions and 
Woods-Saxon wavefunctions. 

IV. RESULTS FOR {ix-,u) TRANSITION RATES 

In Table I we list the measured capture rates and calculated capture rates for the relevant 
(/X, u) transitions to the low-lying "^^K states. We give the results of the WBMB and SAS 
calculations and the Measday and Igo-Kemenes experiments. Note that the experimental 
technique of 7-ray detection in Refs. [16,27] was insensitive to the ground state transition 
and the 30 keV state transition and that Igo-Kemenes et al. only quoted results for the 800, 
1959 and 2047 keV states. Also note that the hw WBMB calculation implies a vanishing 
rate to the positive-parity states at 1632, 1959, 2260 and 2290 keV and the SAS calculation 
gives no prediction for the (^3/2^3/2 states at 2047, 2070, 2104 and 2626 keV. 

Overall - for the ten measured rates to the low-lying states - we claim fair agree- 
ment between experimental data and model calculations from the perspective to the global 
distribution of the capture rates, i. e. the presence of strong, moderate and weak transitions. 
Clearly though - both in the comparison of experiment versus theory and model versus 
model - some discrepancies are obvious. We give below the detailed assessment of the cap- 
ture rates for the various allowed, first-forbidden and higher-forbidden transitions In partic- 
ular, we highhght the (0+, 0)^(0", 2626) axial charge transition, the (0+, 0)^(2", 800) and 
(0+, 0)^(2-, 2047) unique first-forbidden transitions, and the (0+, 0)^(1+, 2290) allowed 
Gamow- Teller transition. 

A. 0+ 0~ first-forbidden transition 

The capture rate for 0+^0^ transitions is determined by two nuclear matrix elements; 
the axial charge matrix element Mqct • V, which originates from the time component of the 
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axial current, and the £ — 1 retarded Gamow Teller matrix element Mqi • a, which originates 
from the space component of the axial current. The interest in MqCt- V stems from its strong 
enhancement by pion exchange currents while the interest in Mqi -a term stems from the large 
contribution of the induced couphng Qp. Recent studies of the ^^O(0+, 0)— >^^N(0~, 120) axial 
charge transition within large-basis shell model calculations have yielded values ofgp — 6 — 9 
[40] and gp = 7.5 ±0.5 [41] that are consistent with the theoretical prediction gp = 8.23 from 
chiral symmetry arguments [9-11]. 

Our calculated rate for the (O"*", 0)— i^(0~, 2626) transition as a function of the couphng 
gp is shown in Fig. 2. The calculations were performed with the hcu WBMB model and 
the (0+2) hu; WBMB model. Fig. 2 also shows the results for (i) only the axial charge 
matrix element Mqu • V, (ii) only the £ — 1 retarded GT matrix element Mqi • a, and (iii) 
both nuclear matrix elements. Note that the curves (i), (ii) and (iii) indicate the similar 
magnitude and constructive interference of the two contributions {i.e. Mqi -a and Mqct- V) to 
the (0+, 0)^(0-, 2626) capture rate. Using the experimental value of A = (15 ± 2) x 10^ s"^ 
for the (0+, 0)^(0", 2626) rate the huu WBMB calculation yields a value of gp = U.S^H 
and the (0+2) hu; WBMB calculation yields a value of gp — lO.Sl^ig. The former value is 
significantly larger than, and the latter value is marginally consistent with, the theoretical 
prediction gp = 8.23 [9-11]. 

Note our calculations have omitted the effects of the exchange current renormalization 
of the axial charge matrix element. Such effects are expected to enhance the matrix element 
Mo(T ■ V by typically a factor of two (sec Rcf. [42] for details). A ~50% enhancement of 
the Moa ■ V matrix element in the (0+, 0)^(0~, 2626) transition would further increase the 
calculated capture rate and thereby further increase the induced coupling gp. 

A crucial issue in the determination of the coupling gp from the (0^,0)^(0^,2626) 
rate is the model uncertainties in the ^°Ca g.s. wavefunction. As the multi-particle, multi- 
hole component of the ^'^Ca ground state is increased then the muon capture rate for the 
(O"*", 0)— >(0~, 2626) transition is decreased, this explaining for the lower capture rate and 
derived g'p-value in the (0+2) tiw WBMB calculation and the higher capture rate and de- 
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rived g'p- value in the hu WBMB calculation. As discussed later in Sec. IV C, we find 
that neither the fkv WBMB calculation nor the (0+2) huj WBMB calculation give a sat- 
isfactory description of the multi-particle, multi-hole admixtures in the *^°Ca ground state. 
Unfortunately, we believe this result forestalls a firm conclusion on the in-medium value 
of the induced pseudoscalar coupling Qp from the (O"*", 0)— i^(0~, 2626) capture rate (rather 
we suspect that the ^°Ca g.s. wavefunction uncertainties are the most likely cause of the 
unexpectedly large value of the coupling Qp derived from the (O"*", 0)— >(0~, 2626) transition 
rate) . 

We note the effects of np-nh admixtures in first forbidden beta-decay were considered 
by Warburton et al. [23]. These authors concluded that the multi-particle multi-hole com- 
ponents of the nuclear wavef unctions, while yielding relatively large corrections for AJ'^ — 
2" first forbidden transitions, would yield relatively small corrections for AJ'^ = 0~, 1~ first 
forbidden transitions. Specifically, they found an increase of the Mqu • V axial charge matrix 
element and a decrease in the Mqi • a ^-retarded GT matrix element of roughly 10%, a small 
correction that would largely cancel in the (0"*", 0)— >(0~, 2626) capture rate. Herein, the ap- 
parent differences between the hu WBMB calculation and (0-1-2) hu WBMB calculation, 
and different closed-shell occupancies in the various model calculations, is indicative of a 
significantly larger model uncertainty. 



B. 0+ — >^ 2 unique first-forbidden transitions 



The experimental data also include two unique first forbidden O"*" — > 2^ transitions to 
low-lying states - the (2~, 800) member of the d^^2~h/2 quadruplet and the (2~, 2047) 
member of the c?3/2"P3/2 quadruplet. The model calculations indicate both transitions are 
governed by the weak axial coupling ga and the i = 1 matrix element M2i-a, with the 
(0+, 0)^(2~, 800) transition being dominantly (^3/2— >/7/2 in character and the (0+, 0)^(2~, 
2047) transition being dominantly (^3/2— >P3/2 in character. 
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Concerning the (0+, 0)— >(2~, 800) transition both model calculations and experimental 
data are in qualitative agreement on the very large capture rate. However, in detail the 
calculated rates are larger than the measured rates by factors that range from ~2 for the 
(0+2) hu; WBMB model to ~3 for the two other models. Concerning the (0+, 0)^(2", 2047) 
transition both model calculations and experimental data are in qualitative agreement on the 
rather moderate capture rate. Unfortunately the experimental results are themselves only 
marginally consistent, with the earlier result of Igo-Kemenes et al. being entirely consistent 
with the calculated rates while the later result of Measday and Stocki being somewhat larger 
than the calculated rates. 

For the first forbidden (0+,0) — >(2", 800) transition the analog transitions in electron 
scattering ^0Ca(e,e')^°Ca(8.43 MeV) [17-19], proton scattering 40Ca(p,p')^°Ca(8.43 MeV) 
[43] and charge exchange '^°Ca(p,n)^°Sc(0.77 MeV) [20,21] are all well known. In Table II 
we compare the 0"^— >2~ reduced transition probability B(M2) obtained from the ^'^Ca(e,e') 
data [19], and the 0+— >2~ capture rate A obtained from the ^°Ca(//, z/) data [16], with the 
corresponding predictions of our model calculations and a pure ^3/2—^/7/2 single-particle 
estimate. The ^3/2—^/7/2 single-particle estimates are seen to exceed the measured value of 
B(M2) by a factor 6.3 and the measured value of A by a factor 4.6. By comparison the model 
calculations are seen to exceed the measured values of B(M2) by factors of 1.5-3.1 and the 
measured values of A by factors of 1.6-2.8. The similar quenching for the 0"*"— >2~ transition 
in the muon capture data and the electron scattering data is not surprising, both processes 
being dominated by the matrix element of the spin-dipolc (M2i-(t) operator. Comparable 
quenching is also reported for the analog transitions in the ^°Ca(p,p') data and the ^°Ca(p,n) 
data. 

This quenching or hindrance of the measured rate over the calculated rate - as observed 
for the {iJtr,v) capture rate and the (e, e') reduced transition probabihty to the (2~,800) 
state - is well documented for unique first-forbidden /3-decay in mass ^4 ~ 40 nuclei (for 
example see Warburton et al. [23]). As documented in detail by Towner et al. [44], such 
quenching originates from the repulsive character of the T — 1 particle-hole interaction and 
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its coherent effects on the particle-hole admixtures in the initial-final state wavef unctions. 
In particular, in Ref. [23] the authors report a typical hindrance of approximately 3.7 for 
the known unique first-forbidden /3-decay transitions in the ^4 ~ 40 mass region in a global 
comparison to their WBMB calculations. 

By comparison we observe no hindrance of the measured capture rate compared to the 
calculated capture rate for the {0^ ,0)—>-{2~ , 2047) transition. Interesting, the (0"'",0)— >(2~, 
2047) transition is dominantly (^3/2— >P3/2 in character whereas the (0'^,0)— >(2~, 800) transi- 
tion is dominantly d^/2^f7/2 character (the unique first- forbidden beta decays in the A ~ 40 
mass region are also dz/2^h/2 transitions). As described by Towner et al. [44], the hin- 
drances associated with T — 1 particle-hole interactions are dependent on the orbital angular 
momenta of the particular transition, and might differ between the c?— >/ angular momenta 
of the (0+,0)— >(2~, 800) transition and the d^p angular momenta of the (0"'",0)— >(2~, 2047) 
transition (of course it would be premature to draw any such conclusion from a data-set of 
two 0+— >2~ transitions). 

C. 0+ — > 1+ allowed Gamow- Teller transition 

The leading contribution to 0+— >1+ transitions originates from the weak axial coupling 
constant Qa and the spin-flip matrix element MiQ-a. Since the MiQ-a matrix element vanishes 
in a simple closed-shell description of the ^°Ca ground state such 0+— >1+ transitions are 
especially sensitive to the np-nh admixtures in the ^'^Ca g.s. wavefunction. 

As discussed in Sec. II C the speciflc np-nh admixtures in the "^"Ca g.s. are strikingly 
different in the SAS model, huj WBMB model, and (0+2) huj WBMB model. In the 
SAS calculation the (0+, 0) — > (l"*", 2290) transition is dominated by a large c?3/2^c?3/2 one- 
body transition density that originates from the large d^^^-fy^ admixture in the "^^Ca g.s. 
wavefunction. In the (0+2) hw WBMB calculation the (0+,0) (l+,2290) transition is 
comprised of many small one-body transition densities that originate from the fragmented 
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np-nh admixtures in the ^'^Ca g.s. wavef unction. Consequently, the calculated (0+, 0) — > 
(l+,2290) rate is relatively large in the SAS model (19.2 x 10^ s~^) and relatively small in 
the (0+2) huj WBMB model (0.08 x 10^ s"^). The measured rate (13±5) x 10^ s'^ is slightly 
smaller than the SAS prediction but much larger than the (0+2) hu; WBMB prediction, thus 
suggesting a better treatment of the np-nh admixtures in the ^°Ca g.s. wavefunction by the 
SAS model. 

This conclusion is supported by the inelastic electron scattering data for the 
^°Ca(e,e')'^°Ca(l+, 9.87 MeV) analog transition. This transition was studied by Petraitis 
et al. [19] and yielded a reduced transition probability of B(M1)= 0.32 ± 0.09 ii% to be 
compared with a SAS prediction of B(M1)= 0.86 //^ and a (0+2) hu WBMB prediction 
of B(M1)~ 0.01 The moderate over-prediction of B(M1) by the SAS model and gross 
under-prediction of B(M1) by the WBMB model is quite similar to the corresponding com- 
parisons between the calculated rates and the measured rates in the ^l capture reaction. 
The similar scaling between experimental values and model values in fx capture and (e,e') 
scattering is not surprising as both processes are dominated by the spin-flip MiQ-a matrix 
element. 

As a final demonstration of the high sensitivity of the (0+, 0)— *-(l+, 0) transition to the 
np-nh components in the /I = 40 wavefunctions we list in Table III the (0+, 0)— >(1"'", 2290) 
capture rate versus the maximum occupancy of the /7/2 orbital in the SAS model. It clearly 
shows the strong correlation between the 0+— >1+ transition rate and the np-nh admixtures 
in the '^'^Ca, 0+ ground state and the '^'^K, 1+ excited state. For example, by increasing 
the SAS model space from 0-2 huj to 0-8 Tiuj one increases the O+^l"*" capture rate from 
4.7 X 10^ s-^ to 19.2 X 10^ s^^. 

D. Other transition 

The remaining transitions comprise the (0+, 0)— >(0+, 1632) non-analog Fermi transition. 
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the (0"*", 0)— 2104) mixed first-forbidden transition, and four £ > 1 forbidden transitions 
that range from second-forbidden to fifth-forbidden. 

The leading contribution to (0"*", 0)— >(0"'", 1632) transition originates from the weak vector 
couphng constant and the nuclear matrix element Mq. As discussed earlier the transition 
rates for low-lying positive-parity states are highly sensitive to the multi-particle, multi-hole 
admixtures in the ''^Ca ground state. In the SAS calculation the transition involves the 
destructive interference of d^/2-^dz/i and /7/2— *-/7/2 s.p. transitions and yields a capture 
rate 1.2 x 10^ s~^. In the (0+2) Uoj WBMB model the transition involves many small 
contributions from different s.p. transitions that destructively interfere and yield a capture 
rate 0.01 x 10^ s~^. Like the (0+, 0)— >(1+, 2290) transition the large differences in calculated 
rates are a consequence of the differences in the multi-particle, multi-hole admixtures in 
the "^"Ca g.s. model wavef unctions. The measured rate of (13 ± 10) x 10^ has large 
uncertainties and is consistent with both calculations. 

The leading contribution to (0+, 0)— >(1~, 2104) transition originates from the weak axial 
coupling constant Qa and the nuclear matrix element Mn-a. The transition to this d^/^2~P3/2 
particle-hole state is dominated by the (^3/2— >P3/2 single particle transition. The measured 
rate of (18 ± 5) x 10^ s'^, ku WBMB rate of 14.7 x 10^ 5-1, and (0+2) huj WBMB rate 
of 9.0 X 10^ s~^, are consistent within 1-2 standard deviations. 

The various model calculations for the four £ > 1 forbidden transitions all imply small 
capture rates. While the measured rate for the fifth-forbidden 0"^— >(5~, 891) transition 
and the third-forbidden 0+— >(3+, 2260) transition are reasonably consistent with model 
predictions, the third forbidden 0+^(3", 2070) transition and second forbidden 0+^(2+, 
1959) transition have unexpectedly large measured rates for ^ = 2,3 forbidden transitions. 
We speculate - as mentioned by Measday and Stocki [16] - that their experimental capture 
rates may include unidentified cascade feeding from higher-lying '^°K levels. 
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V. RECOIL LONGITUDINAL ALIGNMENT 



Both 0+— ^l"^ and 0+— >2~ transitions involve contributions from the induced pseudoscalar 
couphng constant Qp. While the contribution of Qp to their capture rates is generally quite 
small the contribution of Qp to their longitudinal alignments are frequently rather large.^ 
Given the observation of 0+— >1+ and 0+— >2~ transitions in ^°Ca(/x,i/) capture, and deter- 
minations of the analogous transition probabilities in ^°Ca(e, e') scattering, we thought it 
interesting to consider the possibility of determining Qp by measuring the alignment 
following ^°Ca capture. Such measurements were previously performed for muon capture on 
i^N [45], 28Si [46] and ^^Cl [47]. 

Figs. 3 and 4 give results for the longitudinal alignment 02 versus the induced coupling 
Qp for the 0+ (2-,800) transition and the 0+ (l+,2290) transition using the SAS 
model (we adopt the SAS model rather than the WBMB model as it reasonably reproduces 
both the capture rates). Shown are results using (i) all the contributing nuclear matrix 
elements and (ii) only the leading matrix elements (either Miq-o" or M2VO'), In addition, we 
plot the results with maximum /7/2 occupancies of either four nucleons or eight nucleons in 
order to gauge the model dependences. Concerning the 0^—>-2~ transition the calculation 
shows a relatively low sensitivity to both the coupling Qp and the nuclear wavef unctions, 
the wavefunction insensitivity reflecting the dominance of the f7/2^d^/2^^ transition and 
the M2i-(T matrix element. Concerning the 0+— >1+ transition the calculation shows both 
greater sensitivity to the coupling Qp and greater sensitivity to the nuclear wavef unctions, 
the wavefunction sensitivity reflecting the interference between ^3/2^1^3/2 and /7/2^/7/2 

^The recoil orientation about the neutrino momentum is termed the longitudinal orientation and 
the recoil orientation about the muon spin is termed the average orientation. For J > 1 recoils the 
orientation includes both a recoil polarization (rank-one orientation) and a recoil alignment (rank- 
two orientation). We consider the longitudinal alignment 02 which is experimentally accessible in 
the 7-ray experiments. 
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single-particle transitions and Mio-a and Micr-V matrix elements. 

Could the longitudinal alignments in 0+ (2~, 800) transitions or 0+ — > (1+, 2290) tran- 
sitions be measured? In certain cases the alignment imparts a directional correlation between 
the recoil nucleus and the subsequent de-excitation 7-ray [49] , and is thereby measurable by 
the 7-ray Doppler lineshape. However, such measurements require both a suitable 7-decay 
spin sequence and a short 7-decay lifetime. Specifically, the 7-recoil directional correlation 
W{9) is given by 

W{0)o^l + a2B2iP2{cos0) (1) 

where P2(cos^) is the Legendre polynomial, 9 the angle between the recoil direction and 
the 7-ray direction, and B21 the 7-radiation coefficient (see Ref. [48] for details). For the 
0+^1+ transition the dominant Ml decay (l+,2290)^(0+,1644) [50] yields B21 = l/\/2 
which implies a comparatively high sensitivity to 02. For the 0"^— >2~ transition the dominant 
Ml decay (2~,800)— i^(3~,30) [50] yields B21 — I/VTO which implies a comparatively low 
sensitivity to 02. In addition, the lifetime for the gamma de-excitation must either be 
shorter or comparable to the slowing-down time of the recoil ion. Here the 0"'"^1"'" lifetime 
of 83 fs and the 0^^2~ lifetime of 280 fs are comparable to the recoil slowing-down time of 
~250 fs [51]. In short, we suspect such measurements are possible but challenging. 

VI. SUMMARY 

In summary we report the calculation of partial rates and recoil orientations in muon 
capture on calcium-40. The calculations were performed with well-established, well-tested 
nuclear models: the WBMB interaction and sd-pf model space of Warburton et al. [30] and 
the SAS interaction and (1^/2-^7/2 model space of Sakakura et al. [34]. Moreover, the WBMB 
calculations were performed with both a simple closed-shell ^'^Ca ground state and a mixed 
(0-1-2) fiu; ^^Ca, ground state. Taken together the calculations were capable of reproducing 
the important features of low-lying negative and positive parity ^'^K levels. 

Overall we observed fair agreement between measured capture rates and calculated cap- 
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ture rates for the low-lying levels, i.e. the general distribution of muon capture rates to 
low-lying levels being reasonably consistent between the model calculations and the ex- 
perimental results. We note however two striking exceptions - the third forbidden 0"*"— >(3~, 
2070) transition and second forbidden 0+— >(2+, 1959) transition - for which the measured 
rates exceeded the calculated rates by very large factors. 

Most importantly we emphasized the (0"*", 0)— >(0~, 2626) axial charge transition and its 
high sensitivity to the induced pseudoscalar coupling gp of the proton's weak interaction. Us- 
ing the measured capture rate of Measday and Stocki [16] our (0+2) ha; WBMB calculation 
yielded gp = IQ.Sllfl 

huj WBMB calculation yielded gp = 14.3li:6. Unfortunately, 
we concluded that neither the simpler ku; WBMB calculation nor the richer (0+2) hu 
WBMB calculation were capable of satisfactorily describing the multi-particle, multi-hole 
admixtures in the "^°Ca ground state, thus forestalling a firm conclusion on the in-medium 
value of the coupling constant gp from the (0+, 0)— >(0~, 2626) transition rate. 

The above deficiencies in the model descriptions of the ^'^Ca ground state were highlighted 
by our discussion of the (0"'",0)— >(!"'", 2290) allowed Gamow- Teller transition. Such 0"*"— ^l"*" 
transitions on doubly-magic ^°Ca are especially sensitive to the np-nh admixtures in the 
^°Ca ground state as the MiQ-a spin- flip matrix element vanishes for a simple closed-shell 
wavefunction. Unfortunately, we found large model-to-model variations in the different 
calculations of the (0+, 0)^(1+, 2290) capture rate, and most worrisome for the analysis of 
the (0+, 0)— >(0~, 2626) transition, the (0+2) hu; WBMB calculation grossly under-estimates 
the (0+,0)^(l+, 2290) transition rate. 

In addition, we discussed two O"*"— >2~ unique first forbidden transitions to levels 
at 800 and 2047 keV. In nuclear beta-decay such first-forbidden transitions have been ex- 
tensively studied in the context of their hindrance via the coherent effects of the repul- 
sive T—1 particle-hole interaction. Intriguingly, we found a substantial hindrance of the 
(0"*", 0)^(2", 800) transition that has dominant (^3/2— >/7/2 character but a neghgible hin- 
drance of the (0^, 0)— >-(2~, 2047) transition that has dominant c?3/2— *P3/2 character. The 
latter case of a o?3/2— ^P3/2 transition is generally inaccessible via the /3-decay studies in the 
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A AO mass region. 

A number of improvements could be made on our work. The calculations were conducted 
using limited model spaces, harmonic oscillator nuclear wavefunctions, and uniform muonic 
wavef unctions. In particular, the multi-particle, multi-hole structure of ^°Ca ground state 
was a major difficulty for the model calculations, it limiting our conclusions for the induced 
coupling gp in the (0"*", 0)— >(0~, 2626) transition. Given the sensitivity of the 0"*"— >0~ axial 
charge transition to the induced pseudoscalar coupling, the interest in the 0"^— >2^ unique 
first forbidden transitions due to core polarization effects, and the interest in the 0"*"— 
allowed Gamow- Teller transitions due to ground state correlations, we strongly encourage 
further theoretical work on ^°Ca muon capture. In addition, new experimental efforts on 
exclusive muon capture in other ^4 ~ 40 nuclei would help to extend the experimental 
data-set and benefit any model studies. 

We wish to thank both Prof. David Measday and Dr. Trevor Stocki for valuable discus- 
sions and gratefully acknowledge the National Science Foundation (USA) for their financial 
support. 
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TABLES 

TABLE I. Comparison of measured rates and calculated rates for muon capture to low-lying 

levels. The first two columns list the spin-parities and excitation energies of the ^'^K states and 
the final column lists the dominant particle-hole configurations of the ^'^K states. The experimental 
rates are taken from the experimental work of Measday and Stocki (column three) and Igo-Kemenes 
et al (column four) and the calculated rates were obtained with the SAS, Tiw WBMB, and 
(0-1-2) huj WBMB models (see text for details). 
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TABLE II. Comparison of measured and calculated muon capture rates A, and measured and 
calculated reduced transition probabilities B(M2), for the (0+,0)— >^(2~, 800) unique first forbidden 
transition. Column three lists the ^3/2^/7/2 single-particle estimate and columns four, five and 
six list the results of the SAS, Hu; WBMB and (0+2) Hu; WBMB models (see text for details). 



observable measured ^3/2^/7/2 SAS hu WBMB (0+2) hu WBMB 

value s.p. estimate model model model 



B(M2) (jUiv^ fm2) 235+20 1495 778 692 375 

A (x 10=^8-^) 127+13 600 339 351 200 
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TABLE III. The calculated muon capture rate for the (0+,0)^(l+,2290 keV) allowed 



2p-2h and 4p-4h admixtures in the Ca and K wavefunctions. 



maximum ^°Ca 2p-2h ^°Ca 4p-4h (%) 2p-2h 4p-4h rate 

/7/2 occupancy admixture (%) admixture (%) admixture (%) admixture (%) (xlO'^ s^^) 



2 27.0 100 4.7 

4 35.7 8.9 71.3 28.7 11.0 

6 35.5 13.6 56.5 36.0 16.8 

8 35.0 14.5 52.5 36.3 19.2 
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FIGURES 
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(^3/2)"' (P3/2)' 2.047 2- 
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FIG. 1. Summary of assignments of model states to experimental states for the low- lying ^'^K 
levels. Shown are the spin-parities, excitation energies ,and dominant particle-hole configurations 
for the low-lying ^'^K states. 
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FIG. 2. Calculated rate for the (0+, 0)— >^(0~, 2626) transition as a function of the coupling 
constant gp from the (0+2) hco WBMB calculation (upper panel) and the hco WBMB calculation 
(lower panel). The open circles correspond to the axial charge matrix element only, the filled 
triangles correspond to the £ = 1 retarded GT matrix element only, and the filled circles correspond 
to the full calculation. The measured rate A = (15 it 2) x 10^ s~^ is shown by the horizontal lines. 
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FIG. 3. Calculation of the longitudinal alignment 02 for the 0+— >-(l+,2290) transition versus 
the coupling gp. The solid circles indicate the SAS calculation permitting up to eight particles in 
the /7/2 orbital and the open circles indicate the SAS calculation permitting up to four particles in 
the /7/2 orbital. The calculation employing only the M\Q-a matrix element is denoted by triangles. 
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FIG. 4. Calculation of the longitudinal alignment 02 for the 0+— >-(2 ,800) transition versus the 



coupling Qp. The solid circles indicate the SAS calculation permitting up to eight particles in the 



/7/2 orbital and the open circles indicate the SAS calculation permitting up to four particles in the 
/7/2 orbital. The calculation employing only the M2i-a matrix element is denoted by triangles. 
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